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, Recent HBT results from the CERES experiment at the SPS are re- 

p viewed. Emphasis is put on the centrahty and beam energy dependence, 

<^ ' and the results are put into perspective with results at lower and higher 

beam energies. The rather weak beam energy dependence of the HBT 
radii may be understood in terms of a transition from baryon to pion 
dominated freeze-out. The observed short lifetimes and emission dura- 
. tions are presently in contradiction to results from model calculations. 

1 Introduction 

^ . 

O ■ While single particle momentum distributions give only indirect insight to the 

^ ! space-time evolution of the system created in high energy heavy ion collisions, 

^ \ the lifetime and the spatial extent of the pion source as well as the existence 

Q^l of collective velocity fields at the time of thermal freezeout can be probed by 

Q_|' the study of Bose-Einstein momentum correlations of identical pions via HBT 

^ ■ interferometry The width of the correlation peak at vanishing relative mo- 

menta reflects the so-called length of homogeneity of the pion emitting source. 
Only in static sources can the length of homogeneity, in the following also 
^ ' called 'source radius', be interpreted as the true geometrical size of the sys- 

' tem. In a dynamic system, the occurence of space-momentum correlations of 

the emitted particles due to collective expansion generally leads to a reduction 
of the observed source radii, depending on the strength of the expansion and 



the thermal velocity of the pions ^jTfjmt at thermal freeze-out p|. A differen- 
tial analysis of the HBT correlations in bins of the pair transverse momentum 
^t=||p<^i ^Vt,'i\ thus provides valuable information about the properties of the 
collective expansion of the system P|. 

To obtain most detailed information about the space-time evolution the 
three-momentum difference vector q of two like-sign pions is decomposed into 
components, g=(q'iong, fede, Q'out)- Following Bertsch and Pratt 0], giong is the 
momentum difference along the beam direction, calculated in the longitudinal 



rest frame (LCMS) of the pair, gout is parallel to the pair transverse momentum 
kt and gside is perpendicular to giong and gout- 

The correlation function is defined as the ratio C2{q)=A2{q) / B2{q) where 
the 'signal' A2{q) is the probability to find a pair with momentum difference g 
in a given event and the 'background' -82(1?) is the corresponding mixed-event 
distribution. 

The normalized correlation functions are fit by a Gaussian: 
C2(g) = 1 + >^exp{-Rl^^^ql^^- Rl^^ql^^ 

~-^out5'out ~ 2-Routlong5'outQ'long); (1) 

with -Riong, -Rside, -Rout being the Gaussian source radii and A the correla- 
tion strength. The cross-term -Routiong appears as a consequence of space-time 
correlations in non-boost-invariant systems. 



2 Results from SPS 



Pion HBT interferometry analyses have been performed by a number of ex- 
periments at SPS H, 1^, H- A systematic study of the centrality dependence 
of HBT radii around midrapidity from Pb+Au collisions at all three presently 
available SPS energies, 40, 80, and 158 AGeV, was recently presented by the 
CERES collaboration [^]. Details of their analysis can be found in jlOl . 
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Fig. 1: i?iong vs. kt The lines represent fits to the data (see text). 

The fcrdependence of the longitudinal source parameter -Riong is shown in 
Fig. 1 for the three different beam energies and as function of the centrality 
of the coUision. Similar to previous studies, a strong decrease of -Riong with 
kt is observed, characteristic for a strong collective longitudinal expansion, 



where the length of homogeneity is entirely saturated by the thermal length 
scale ^ 



'Tf/mt. For the case of longitudinal boost-invariance, -Riong can be 
connected to the lifetime of the system, the time elapsed between the onset 
of expansion and kinetic freeze-out by -Riong=Ts(7//mf) 2 

The slight but systematic overall increase of -Riong with centrality and beam 
energie is reflected in a correspondingly longer lifetime when applying the above 
expression to the data and assuming Tj=120 MeV (see fits in Fig. 1). The 
observed lifetimes at the SPS are 6-8 fm/c. 
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Fig. 2: ii'.idc vs. kt 0. 
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The transverse source parameter -Rside is closely related to the true geomet- 
rical transverse extension of the system at freeze-out. Consequently, a weaker 
/ct-dependence is observed for -Rside (Fig. 2), indicating the presence of radial 
flow. A slight but systematic increase of -Rside with centrality is measured at all 



beam energies, as expected from simple collision geometry. The beam energy 
dependence is surprising: largest -Rside and the strongest /cj- dependence are ob- 
served at the lowest energy. This led previously to the conclusion that radial 



flow may reach a maximum at the lower SPS energies ITJ]. We will discuss 
below that the extraction of the radial flow velocity from the fcf-dependence of 
i?side at the lower SPS energies might be questionable. At 158 AGeV, a radial 
flow velocity of about 0.5c was obtained from the analysis of single particle 



mrspectra and -Rside (^t) R |T2 



Many authors have discussed that a strong first order phase transition with 
a large latent heat can lead to a retardation of hadronization during the mixed 



phase and consequently to a long duration of pion emission |T^, |T4|, [15 



J' 



T6|. This should be observable via a strong increase of the outward source 



parameter i?out with respect to -Rside- However, no long-lived source has been 
observed so far. At SPS, -Rout (Fig. 3) is similar to -Rside, indicating a short 
pion emission duration and a sudden freeze-out. 



3 Beam energy dependence 

A large amount of pion interferometry data have been published by exper- 
iments at AGS, SPS, and RHIC. This allows for a systematic study of the 
source parameters over a wide range of beam energies. It has been argued 
that the HBT null effect [|17|, the absence of any beam energy dependence, 
in particular when going to RHIC, is the actual surprising result from HBT 
interferometry. 

In Fig. 4 are shown the fcf-dependences of -Riong, -Rside, and -Rout in central 



Pb(Aa)+Pb(Au) collisions at different beam energies |]T8|, |T0|, |T9[. Indeed, no 
dramatic variation of any of the source parameters can be observed. However, 
a closer inspection reveals some interesting features. -Riong is approximately 
constant from AGS to the lower SPS energies, but develops a significant in- 
crease within the SPS regime and towards RHIC, indicating a smooth increase 
of the lifetime. 

Most interesting is the behaviour of -Rside: It is gradually decreasing at 
small kt up to top SPS energy, connected with a continuous flattening of the 
/ct-dependence. At RHIC, -Rside is again larger than at SPS while the shape is 
not yet well measured, at least the 7r~7r~ sample looks rather flat. Naivly, the 
flattening indicates a decrease of the radial flow velocity as function of beam 
energy, in contradiction to the present understanding ||20[ . 



-Rout shows a rather weak energy dependence, possibly indicating a shallow 
minimum at the lower SPS energy. 
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Fig. 4: /ct-dependence of HBT radii around midrapidity in central Pb(Au)+Pb(Au) collisions 
at different beam energies. The data are from |lO[ 0. 



A straight-forward investigation of the freeze-out properties can be per- 
formed by relating the measured source parameters to an effective freeze-out 
volume, V/=27r_Riong-Rsidc- Assuming freeze-out at constant density we 



expect Vf to scale linearly with the charged particle multiplicity. Fig. 5 shows 
Vf as function of the number of participants at 40, 80, and 158 AGeV. A lin- 
ear scaling with A^part is indeed observed at all three energies, consistent with 
the assumption of a constant freeze-out density, since the number of charged 



particles was found to scale close to linear with Apart at SPS |^ . 
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Fig. 5: Vf as function of A'^part at different SPS energies [[fO[ 



"^000 



cr!l500 



DC 

(=; 

1000 



500 



<k,>=0.16 GeV/c 



T E895 

• CERES 

* STAR nn 





1 10 lO'' 

Vs(GeV) 

Fig. 6: Vf in central Pb(Au)+Pb(Au) collisions as function of ^/s (see text). E895 points 
are at < h >=0.15 GeV/c @, STAR is at < h >=0.17 GeV/c |l9|, CERES points are 
the average of their bins at < h >=0.12 and < kt >=0.2 CeV/c ]1C|]. 



But the beam energy dependence is surprising: There is no clear hierarchy 
visible in Fig. 5 as expected from the increase of multiplicity by about 50% 
between 40 AGeV and 158 AGeV; the smallest Vf are observed at 80 AGeV. 
Obviously, the freeze-out volume scales with multiplicity as long as multiplicity 



is controlled via centrality, but it does not scale accordingly as multiplicity 
changes with beam energy. 

The comparison of Vj at different beam energies from AGS to RHIC sheds 
some light on this: The freeze-out volume Vf is gradually decreasing within 
the AGS energy range, reaches a minimum at SPS and then increases towards 
RHIC, as demonstrated in Fig. 6. 

Clearly, a simple relation between multiplicity and freeze-out volume can- 
not hold. On the other hand, it is plausible to assume that at low energies 
pions interact mainly with nucleons, while at high energies pion-pion scatter- 
ing dominates. Is the transition from nucleon to pion dominated freeze-out 
characterized by a minimum in Vf at SPS? 
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Fig. 7: Left: Midrapidity yields of different particles vs. y/s. Upper right: cross section 
weighted sum of particle yields (see text). Lower right: ratio Vf/Ncff vs. y/s. 

In Fig. 7 (left) the midrapidity density of pions |2^, ^ ^ and pro- 
tons |23, 0, [11], ^, |3g, |31|] in central Pb(Au)+Pb(Au) collisions is shown as 
function of ^/s. The pion yield increases monotonically with beam energy. 
The total proton (p+p) yield at midrapidity drops from AGS to SPS and stays 
approximately constant between SPS and RHIC because the decreasing num- 
ber of net protons is compensated by p-p production. The sum of pions and 
protons, however, is still a monotonic function of ^/s. 

At this point, the different cross sections a-^^ and (T^at have to be considered. 
In Fig. 7 (upper right) the cross section weighted sum of pions and protons 
A^'eff=2 • iVp+p ■ (J^Ar+ 3 ■ Nt^- ■ (TTTTr is showu [|. For the cross sections (7^,^=10 mb 

"'^Note that there are a few simplifications: the protons arc multiplied by 2 to account for 



and (T^Ar=65 mb are assumed. Indeed, the cross section weighted sum of pions 
and nucleons at midrapidity is non- mo notonic and exhibits a minimum at SPS. 
As a consequence, the ratio Vf/Nes, which has the dimension of a length, is 
approximately constant (Fig. 7, lower right). This result suggests Vf/Nes^l fm 
as a universal freeze-out condition. 

In this picture, the relatively weak beam energy dependence of HBT source 
parameters can be understood as an interplay between the decreasing (and 
eventually levelling off) (p+p) yield and the increasing pion multiplicity at 
midrapidity, if their different cross sections with pions are taken into account. 
It is interesting to obtain a detailed understanding of the role of protons for 
pion freeze-out at low energies. The strong momentum dependence of the cTtttv 
cross section may affect the /ct-dependence of -Rside at low beam energies and 
cause the flattening of the fcj-dependence with increasing beam energy, as the 
importance of protons ceases. At the lower SPS energies, where protons are 
still important, the interpretation of the fcj- dependence of -Rside in terms of 
radial flow may be questionable, and possibly breaks down at AGS. 



4 Discussion 

Single particle spectra and azimuthal anisotropies at SPS and RHIC have been 
well reproduced by state-of-the-art hydro+cascade calculations |3^, ^ |3l| . 
However, they fail to reproduce the measured HBT parameters, in particular 
the lifetime of the system is grossly overestimated. 

In these models, the early plasma phase is described by hydrodynamics 
assuming an ideal plasma EOS with a first order phase transition and a mixed 
phase with an adjustable latent heat. At the end of the mixed phase, hadrons 
are created assuming a thermal phase space population superimposed by the 
collective velocity field created by the plasma pressure. After hadronization, 
hadrons are rescattered using a hadronic cascade code. 



In p5| the response of the hydro+cascade calculation to different choices 
of the latent heat is well documented. From this investigation, it becomes 
obvious that the choice of a relatively large latent heat of 0.8 GeV/fm^, needed 
to describe the SPS single particle data, is mainly driven by the observed mass 
dependence of the inverse slope parameters. Smaller latent heat produces 
larger flow velocities before hadronization and therefore too hard spectra, in 
particular for the multi-strange baryons which do not follow the linear scaling 
behaviour observed for vr, K, p On the other hand, a small latent heat 



the neutrons, which is not completely correct for the lowest beam energies. Also the role of 
light nuclei and other produced particles is neglected. 



would lead to an early acceleration of the matter and therefore to a faster 
dilution and shorter lifetime, in particular of the hadronic phase, which is 
suggested by the HBT measurements. 

This leads to a closer inspection of the data as represented in Fig. 8 (left). 
For some of the particle species, the observed inverse slope parameter scales 
linearly with the particle mass, consistent with the assumption of a common 
flow velocity. There are, however, a number of exceptions which do not follow 
this rule (0, S, Q). This was explained by their small cross sections with the 
surrounding matter, hence they are expected to participate less in the collective 
motion. 
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Fig. 8: Left: Inverse slope parameter vs. mass of hadrons at SPS. Right: Inverse slope 
parameter vs. number of constituent quarks. 

On the other hand, it is evident that the data separate into two groups: 
All mesons have similar slopes, and the same is true for baryons. One may 
therefore as well correlate the inverse slope with the number of constituent 
quarks inside the hadron, rather than with the hadron mass (Fig. 8, right). 
Also here a scaling behaviour can be observed, and almost all particles are 
consistent with this trend. Note that the large inverse slope of the observed 
by NA49 is experimentally still under debate, since there is also a mea- 
surement by NA50 [Q, giving a much smaller number. 

The physical picture which arises from this representation is different from 
the common interpretation: collective motion has completely developed before 
hadronization, with constituent quarks being the flowing objects. Hadrons are 
formed by coalescence of constituent quarks, preferentially if they are close 



in phase space, thereby adding their momenta. When hadrons acquire mass, 
momentum is conserved which leaves the momentum spectra unchanged. In 
this picture, the data are consistent with a kinetic freeze-out temperature of 
0.12 GeV, a collective quark flow velocity of ^1/3 and an effective quark 
mass of 0.33 GeV: T=0.12 + ^ ■ 0.33 ■ nqu^rks (dashed line in Fig. 8, right). 
Most of the explosive power of the system would be assigned to the early, 
pre-hadronic phase of the collision, thereby qualitatively explaining the short 
lifetimes observed by HBT. It also suggests that the hadronic mass scale is 
no more relevant in the early phase of the collision, as naively expected for a 
deconfined partonic system. 
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